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A possibility of formation of static dual scalar and pseudoscalar density wave condensates in 
dense quark matter is considered for the Nambu-Jona-Lasinio model in an external magnetic field. 
Within a mean-field approximation, the effective potential of the theory is obtained and its minima 
("^ ' are numerically studied; a phase diagram of the system is constructed. It is shown that the presence 

Cn ' of a magnetic field favors the formation of spatially inhomogeneous condensate configurations at 

low temperatures and arbitrary nonzero values of the chemical potential. 
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I. INTRODUCTION 



^ . At present, one of the most commonly used effective theories of quantum chromodynamics is the Nambu-Jona- 

^ Lasinio (NJL) model |l|,[3l; a local relativistic four-fermion interaction theory. The QCD and NJL Lagrangians possess 
f^ , the same symmetry group, the NJL model is therefore widely exploited in studying the nonperturbative QCD vacuum 
and its properties under various external conditions. Many features of quarks and light mesons can be successfully 
described within the NJL model on the basis of the spontaneous chiral symmetry breaking phenomenon |3|-|5| . 
J> ■ Considering the QCD ground state properties, a number of studies were dedicated to the possibility of formation 

^^ ' of spatially nonuniform phases in dense quark matter. It was first shown [6| that spatially inhomogeneous and 
OO . anisotropic chiral condensation may occur in QCD at asymptotically high values of the chemical potential and large 
CT^ ■ Nc, the ground state spatial structure taking the form of a standing wave. This phenomenon was further discussed 
^^ in literature [3-Q investigating the possibility of such type of symmetry breaking and its competition with color 
t — ' superconductivity under various conditions including intermediate densities of quark matter. The problem has also 
^D [ been examined recently in the context of quarkyonic matter [10]. Along with QCD studies, similar behavior of the 
^^ ' ground state has also been discovered and successfully reproduced in N JL-like effective models 0, llll - [27| , although 
the effect tends to be dependent on the adopted regularization scheme (see, e.g., Ref. [28| for details). 

Spatially nonuniform condensates proposed at the start of the theoretical research on the subject 11] and stud- 
ied extensively later on are known as dual chiral density waves (DCDW, the name introduced in Ref. ^20!]). The 
corresponding configuration can be described as follows: 



(V'V') = A cos qr, 
(ipi'y^T^'ip) = A sin qr, 

where A is the chiral density amplitude, q is a wave vector (which has to be determined dynamically along with A), 
and Ta are the isospin Pauli matrices. Expectation values (V'V') and (ipij^T^ip) are identified with a and 7r° condensates; 
one generally assumes {ipi^^Tiip) = {ipi^^T2ip) = 0, thus charged tt condensates being absent. In general, scalar and 
pseudoscalar condensates are on the chiral circle: (f/'V')^ + (ipi^/^T^p)^ = A^. It is argued (see, e.g., Refs. [13,1201) ^^at 
DCDW may arise between the massive and symmetric phases of the NJL model at low temperatures if the coupling 
constant is sufficiently large. The formation of DCDW along with color superconductivity has also been discussed 
in literature |2l| - [23| . It should be noted however that, although the majority of studies of condensate inhomogeneity 
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focus on wavelike configurations and DCDW in particular, this is mainly for technical reasons. There may exist other 
competing and even more preferable spatially nonuniform ground state configurations like domain walls, see, e.g., 
Ref. [23|! but, in general, they are much harder to deal with. For the same technical reasons, one usually considers 
the limit of vanishing quark current masses (chiral limit), although recently efforts to get rid of this assumption have 
been made [M~M]- 

The chiral condensation phenomenon (with spatially homogeneous condensate configurations) has recently attracted 
great attention in the situation when external gauge fields and, in particular, strong magnetic fields are present J3Qrt4C|. 
In fact these fields are common in the physical circumstances where the phase structure of quark matter is of interest, 
e.g., in compact stars or in heavy-ion collision processes [4l|. The effect of magnetic fields on the spatially nonuniform 
chiral condensation is therefore worth investigation. This is important in the context of a targeted experimental 
search for possible condensate inhomogeneity signals predicted, e.g., in Ref. jl^. Some interesting results related to 
the subject have been obtained in Ref. [44I showing that a stack of 7r° domain walls may emerge in the QCD vacuum 
in a strong magnetic field due to the axial anomaly (and concerning color superconductivity in a magnetic field, see, 
e.g., Refs. [43, |4J] suggesting the formation of magnetic domains in cores of compact stars). 

In this paper, we examine the formation of the condensate configuration defined in Eq. ([T]) in dense quark matter 
in the framework of the NJL model in the presence of an external magnetic field and show that the latter favors 
the emergence of DCDW at low temperatures. Limiting ourselves to the chiral limit, we base our calculations upon 
exact solutions of the Dirac equation and use the proper-time regularization method such that our results agree with 
Ref. f2(| in the zero-field limit. 

II. THE MODEL 

We start from the NJL Lagrangian density for a quark field ip with Nf = 2 fiavors (representing the up- and 
down-quarks) and Nc = 3 colors: 

/:N.7L ^ ^ ^^^^ -mc + fi-f"") V + G [(V^^)2 + (V^i7Vi^)2] , (2) 

where G is the coupling constant, fi is the chemical potential, rric is the quark current mass, the covariant derivative 
D = d + iQA with A being the electromagnetic field, and Q the electric charge matrix acting in the flavor space: 

We take 7^ — —i'-f^'f^'-f^^^ and assume T3 to be diagonal, we use the standard (Dirac) representation of the 7 matrices 
throughout the paper; matrix indices are suppressed in our notation when possible. In what follows, we hold m^ = 
assuming that the appropriate dimensional parameters in our model tend to be much greater than nic — 5 MeV. 
The symmetry of the model is therefore SUl{2) x SU]i{2) and it is reduced to [/t3l(1) x Ur^nil) when an external 
homogeneous magnetic field is present (with the field strength pointing in the z direction). 

Using ansatz ([T]), we obtain the Lagrangian density in the mean-field approximation (we only take into account 
Hartree terms here, see a discussion on this subject in Ref. [20J'): 

2 
£MF ^ ^ ^ijD + fij'^ — m (cos qr -f 17^x3 sin qr) ] V' — -tt; , (3) 

where we have denoted m ~ — 2GA. We assume that the system resides in an external magnetic field, the wave vector 
q being parallel to the field strength H, both vectors oriented along the z axis. Such an assumption is reasonable due 
to the symmetry considerations; possible small deviations of q from the preferred orientation along H are taken into 
account further. 

As it is commonly done when considering model (j3|), we use a field transformation ip — >■ e*"^ '^^^^ ip, ip —>■ ip e^'^ '^abx ^ 
where b^ = (0, b), x^ = {t, r), and b — q/2, to remove the spatial modulation from the resulting Lagrangian density £: 

9 

777, 

C = 'i{;{i'yD + n-f° -m + j^T3jb)ip~ -—. (4) 

4Cz 

It should be noted, however, that special care is needed when performing such operations in the presence of background 
gauge fields. To obtain correct results, one should apply, for example, Fujikawa's method |45i] and its generalizations 
for finite fermion field transformations. Fortunately, the path integral measure VipVip remains invariant in our case 
since the quantity t^'^"'^ F^vFap which arises in Fujikawa's exponent, where F is the electromagnetic field strength 
and e is the antisymmetric tensor, vanishes in the absence of an electric field. 

In what follows, we obtain the thermodynamic potential $7 for the model described by Eq. Q and then study 
numerically the minima of J7 with respect to the order parameters m and b. 



III. ONE-PARTICLE ENERGY SPECTRUM 

For later convenience, let us first consider a simplified model for a charged fcrmion (electron) field having no flavors 
or colors with the Lagrangian density 

L = tp {i'yD — m — 7^7&) V': (5) 

where D = d — ieA, e > 0. The term ipj^jbip in the latter expression describes a Lorentz- and CPT-breaking 
background interaction controlled by the axial four- vector b^. This type of interaction arising within the context of 
the Standard Model Extension [4^ has been a subject of extensive theoretical research in recent years (see, e.g., Refs. 
[47h50| |V In this paper, in order to obtain the one-particle energy spectrum of model (O, we use a technique similar 
to that adopted in Ref. [50]. 

Let b = (0,0,6), H = (0,0, iJ), H > 0; we take the electromagnetic field in the Landau gauge: A^ — (0,A), 
A = {0,Hx,0). The modified Dirac Hamiltonian derived from Eq. ([5]) is as follows: 

Hu^aP + 7V - E3&, (6) 

where P = — iV + eA is the gauge-invariant kinetic momentum, a — 7^7, S^ = \iijk<^-'^ , cr'^'^ — 5 [7'^, 7"]- Since 
[{olP ±)'^ , Hd] = where Pj^ — {Pi,P2,0), {aP±)'^ being an observable with an oscillatorlike spectrum, it is easy to 
prove that the eigenfunctions of Hq have a standard general form (see Chapter IV of Ref. [5l| for details): 
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{eHfl\ ^ = ^/^x+-l=, (7) 



where u„(^) are the orthonormalized Hermite functions [we assume u_i(^) = 0] and {ci} are spin-dependent coeffi- 
cients. The quantum number n = 0, 1, ... is the Landau level, p is the momentum component parallel to the magnetic 
field direction, and q is related to the symmetry center xq of the wavefunction ^ along the x axis: q = — .tq eH . For 
each n > and fixed q and p, we have an eigenvalue problem for a 4 x 4-sized matrix K acting on the vector {ci}, 
where 



K = aip± + a^p -f- 7 rn — S36, p± — \/2eHn. (8) 

The quantity q is absent in Eq. ([5]) thus providing the degeneracy of the energy spectrum with respect to it; this 
phenomenon is related to the freedom in placing the particle's orbit in a magnetic field and is preserved for any gauge 
of A. 

Let us now consider a unitary transformation: K = U^^KU where U = e^^^^ e^^ ^^^ = ^(1 + il]2)(l + «7°2]2); it 
yields 

A' = aiP_L + asp -f 7°TO + 7°S3/i-f^, P± ^ P, P = -p±, jlH = b. (9) 

The matrix K formally corresponds to an electron with an effective vacuum magnetic moment moving in an effective 
external magnetic field. The problem for this case has been studied and solved in Ref. [52] (note that the form of 
the coefficients {ci} is independent of the adopted electromagnetic field gauge). The case n = requires a separate 
treatment though, since K is reduced to a 2 x 2-sized matrix Ko acting on the coefficients {ci}, z = 2,4: 

j^ _ fm + b -p 
^o-y _p _^ + ^ 

The eigenvalue problem for Xq can easily be solved. The final expression for the energy spectrum has the form 
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E^,^, = r\\^^'^'+P' + V +2ei?"' " = 1'2,..., (^0) 



e 



^Jm? +p^ + 6, n = 0, 



where ^ = ±1 is the spin quantum number, e = ±1 is the energy sign (when n > 0). When n = 0, one only has two 
(instead of four for n > 0) energy branches distinguished by the number e and the latter has lost its meaning of the 
energy sign in the presence of & 7^ 0. Spectrum (jlOp has been known in literature j53| but the energy shift of the n = 
level has not been shown explicitly in the paper cited. The specific asymmetry between the particle and antiparticle 



energy spectra is due to the CPT-odd nature of the background interaction present in our modeL The phenomenon 
does not manifest itself for free particles since one can compensate the CPT- induced transformation b — 7> — b by a 
spatial rotation. But this can no longer be done in the presence of a preferred spatial direction which is introduced 
with H in our problem. 

The coefficients {ci} which meet the orthonormalization requirement for the eigenfunctions {^ngpC^} ^.re as follows: 
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where 
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Formula (jlip is valid for all n provided that one assumes C — ^ when n = 0. This has a physical reason since the 
quantity 11 is an eigenvalue of the spin operator ^^{Ps — j^m) which commutes with Hd, and 11 ~ e^rn? + p^ at the 
lowest Landau level. 

Now that we have found the energy spectrum and a system of wavefunctions for Hamiltonian ([S]), we may use the 
perturbation theory to take into account possible small deviations of b from the direction of the magnetic field. Let 
b = (6_L, 0, &), the corresponding correction to -ffo being Y = —ll\h]_. There is no first-order correction to the energy 
due to the rotational symmetry of the system. The second-order correction obtained through the standard procedure 
is as follows: 
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(12) 



where E is given by Eq. ([T 



R± = V2bi_ l±Ce\/l- 



E^ 



and we have used a stroke symbol to denote that a quantity is a function of the quantum number set {n'p'(^'e'} 
instead of {np(e}, the latter being fixed for a given one-particle state. The numbers n' , e' are expressed through n, e 
differently for each term in Eq. P^ while one has p' — p and ^' = — C in the whole expression. The last two terms 
with n' = n— 1 are absent in the case n = due to R- — provided that C = £• 

Despite the emerging energy level degeneracy with respect to C, Eq. (IT^ is valid in the limit 6 — > 0. It is easy to 
notice though that the terms with e' = e suffer from divergence due to a level crossing possible for states with adjacent 
n and opposite Q, thus making the result obtained not applicable in the corresponding region of the parameter space, 



namely, when 46 (to 



P 



(eH)^ [if ( = sgnb the first term in Eq. p^ is divergent and if C = — sgnfo such is the 



third]. To workaround this, one has to modify the method of calculating AE± in that region. Using the perturbation 
theory formalism for two near-degenerate levels with energies E, E' (see, e.g., Ref. [5J| for details), we find new energy 
values E± with a gap induced by the perturbation V: 



E.^i{E,E'.viE-Er.^RR'). ««'^^r^:: ::;:!;: 



(13) 



Assuming E± = E + AE± and taking into account that 



sgn(^ - E')y/{E-E'Y +ARR' ^ E - E' 



2RR' 



E-E' 
when RR' <C \E — E'\, Eqs. (|12p and ([T3| can be combined into one asymptotic formula with the change 

^^ -^l[-E + E' + sgn(i? - E')^iE-E'r+ARR 

applied to the first and the third term in Eq. (IT^ . with RR' being R+R'_ and R'^R- respectively. The factor 
sgn(_E — E') is used in the latter expression to select a proper branch of solution (IT51) since we want to retain the 



meaning of AE± being a small correction to a particular energy level E when \E — E'\ is not vanishing. There is an 
ambiguity in this approach arising when E = E' in either of the terms which have undergone the change, it may be 
fixed with the help of the following convention: 



sgn(£: - E') 



sgn^{E-E'), n' = n+l, 
sgn_{E-E'), n' = n-l, 



where sgn^(O) = ±1. It is easy to see that this ensures the consistency of the formula obtained (no values of AE± 
have been lost when considering the energy spectrum as a whole). The final result reads 



AE^ = 
1 



i (-^ + E' + sgn+(S - E')^{E-E'Y+AR+R'_ 
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n — n+ 1 



C' — C 



(14) 



The spectrum E + AEj_ with E and AEj_ provided with Eqs. ([TU]) and p^ can now be used to evaluate the effective 
action of the model. 



IV. EFFECTIVE POTENTIAL AND REGULARIZATION 



Let us now return to model @. The corresponding one- loop effective action 



d^x 



- I + - InDet (i'jD + ^7" - m + 7^X376) 



(15) 



is decomposed trivially into similar parts calculated separately for each flavor and color; moreover, it can be expressed 
in terms of the effective action V for the model studied in the previous section with an appropriate change in the 
electric charge and the chemical potential included: 



r = I d'^x\ - 
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4G 
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- InDet (i5" + Ai--ffD), 



(16) 



where Hj^ is given in Eq. ([6]). We have used a charge conjugation for the up-quark when deriving the foregoing. 
Since we know the eigenfunctions {^rigpC«} ^^^ ^^^ spectrum {Enp^e} of Ho, the expression for V can be evaluated 
through the standard procedure: 



r' = ^Jl■ln[-{^^Y + {Hr,-^i)'] 



i/*»i:i„[-(p»)V(E-rt=]^^ij, 



-ip t^ 



where we have introduced a characteristic four-volume LfL^LyLz and 



E = E / ^^ / ^^ = E / '^^ eiJL^. 
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In order to obtain the thermodynamic potential fl, we employ Matsubara's technique |55l [: 
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where /? = 1/T is the inverse temperature; the sum over k is easily evaluated. We finally find 
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Separating the effects of nonzero temperature and the vacuum contribution, expression (J17p can be decomposed 
into three terms: 



n' = n'+ n'„ + n'n 



Tj 



where 



^[, 



n[ 



1 eH 

'2(2^ 

1 eH 



dpJ2{\E- 

nC,e 



t^\-\E\). 



2 (27r)2 



(18) 
(19) 
(20) 



n^e 



The vacuum term $!(, is divergent while the terms J7' and $7^ are finite (being zero when /i = and T = respectively). 
The N JL model is known to be sensitive to the choice of a regularization scheme due to the nonrenormalizable nature of 
the four-fermion interaction [5^] (see also a discussion on this subject for the case of a spatially nonuniform condensate 
in Ref. [2^). We here employ the proper-time method 57 1: 
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where A is the regularization parameter; so our results should agree with those obtained in Ref. 20] in the limit iJ — > 0. 
It is easy to see that expression (1^ for fi^ is well defined, and although expression P^ for f]^ seems to be 
convergent due to the internal sum over e, care is needed when evaluating it since it has been obtained as a difference 
of two divergent objects. One can derive an arbitrary value for such an expression rearranging the terms during the 
summation procedure |58|, so an intermediate regularization is needed to get a correct result. Let it be a simple 
cutoff: 






(21) 



where A' is sufficiently large (not necessarily being equal to A). If the relation E\g=^i = — £'|e=_i > holds, the 
cutoff factor can be dropped out (provided that /i < A'): 

J2 i\E Ml - m e{A' - \E\) ^[{\E-fA + E + f,- 2E) 6{K' - E)] |^^^^ 

e 

= [2 (m - E) 0(/i - E) 0{A' - E)] 1^^^^ = [2 (p - E) 0{^i - E)] l^^^ . 

But this is not the case when the symmetry between the particle and antiparticle spectra is broken: E\^^^i ^ —E\^^-i, 
which occurs at the lowest Landau level in our problem [and when taking into account corrections (|14p to the energy 
levels as well]. In general, one has to retain the regularization throughout the calculations or modify the whole 
expression by a finite but nonzero correction. The effect for the case 6j_ = 0, n = can be studied exactly (see the 
Appendix) : 

[dpY,{\E-fA-me{A'-\E\) = fdpJ2i\E-^^\-\E\) + 4^^b. 

•> g A'— >CXD J g 

If one omits the term 4/i6 in the above expression, the resulting potential il turns to be dependent on b when m — 0, 
and this is physically incorrect according to definition ([T]); no observable quantity may depend on the wave vector of 
a condensate wave with a zero amplitude. 
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FIG. 1: Phase diagrams for (a) supercritical and (b) subcritical values of the coupling constant G at zero temperature. All 
quantities are dimensionless. There is a symmetric massless phase A with no chiral condensate and two chirally broken massive 
phases B and C, the latter being a phase with a nonzero matter density p, whereas p = in phase B. Massive phases B and 
C are spatially nonuniform when H > 0. There is also a new phase D with a strong condensate inhomogeneity (retaining the 
presence of DCDW in the H ^- limit). The phase transitions are first order. There is a crossover between phases C and D 
in a magnetic field strong enough; we have plotted the boundary between them with a dotted line in that region. 



PHASE DIAGRAM 



To construct a phase diagram of the system, we have studied the minima of the regularized thermodynamic potential 
fl numerically vi^ith respect to the order parameters m and b for different values of the chemical potential /i and the 
magnetic field strength H . We tried to find a global minimum in the case of several minima present on the fl surface. 
We used spectrum (fTO)) and took into account corrections (fT4)) to study the stability of the results with respect to small 
deviations of b from the direction of the magnetic field (taking A' = 10 A). We only used dimensionless quantities 
throughout the calculations with A being the characteristic energy scale. In what follows, we denote these quantities 
with the same symbols as the original ones, e.g., m stands for to/ A, etc. We performed integration over the quantum 
number n instead of summation when eH <§C 1, thus being able to consider the limit 7J — > with no singularities. 
The estimate of the maximum relative and absolute error was set at the level of 10^"^ and 10^®, respectively. We 
take the values of /i and VeH from up to 0.8; it should be noted that there is no physical sense in considering high 
values of these parameters since 1 is the (dimensionless) regularization constant in our model. The critical value of 
the coupling constant is Gc — 3.27 in our model, li H = and // = 0, spontaneous chiral symmetry breaking only 
occurs when G > Gc', we take this fact as the definition of Gc- 

At present, an exact form of the one-particle energy spectrum in the case of 6_l > {m > 0) is not known, so 
that comprehensive analysis of the problem cannot be made. There is no strict guarantee that there are no global 
minima of the thermodynamic potential somewhere in the region to > 0, 6 > 0, &_l > when 7J > 0, /i > since 
one can construct a dimensionless ratio y/eH/fi and the latter may be related to the ratio b±/b. If this is true, the 
DCDW wave vector orientation would be diverted from the preferred direction of H and the rotational symmetry 
of the system would be completely broken. Nonetheless, it is reasonable to believe that the global minima of $7 are 
reached when b± ^ implying that the rotational symmetry is still preserved. To test this to the extent possible, for 
each minimum found (when H > 0), we studied the behavior of the thermodynamic potential in the region of b± close 
to zero. We actually calculated the second derivative d^il/db^lb^^o numerically, and we made use of the explicit 
energy spectrum corrections (fT4|) during the evaluation of that quantity. The latter turned to be positive everywhere, 
so, in this approximation, no instability of the thermodynamic potential minima with respect to b± has been found. 

The results of numerical analysis in the case of T = and supercritical G = 6 are presented in Figs. [T^, [21 |3l As one 
would expect, we recover the result obtained in Ref. [20] in the limit il — >■ 0, see Fig. [3]Ja); and there is a nontrivial 
behavior of the system when H > 0. The order parameter b related to the DCDW wave vector (6 — q/2) grows either 
smoothly (for the range of the chemical potential /x up to some value) or discontinuously (for higher values of /x) with 
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FIG. 2: The order parameters as functions of the magnetic field strength for various values of the chemical potential at zero 
temperature, and G = 6. All quantities are dimensionless. 



the increase of the magnetic field strength i?, the effect being more vivid for greater /x, see Fig. [2] There is also a gap 
corresponding to a transition from a symmetric phase present in a weak field in dense matter (we assume 6 = when 
TO = although h has no physical meaning in that case and can be set to have an arbitrary value). DCDW is absent if 
/i = 0. Noticeably, the order parameter h grows linearly with the increase of the chemical potential \x (up to a critical 
value where a phase transition occurs), the growth rate being higher in the stronger field, see Fig. |3l DCDW is absent 
when H — except for the range of /i corresponding to a new phase examined in Ref. (20| . This phase (we name it 
phase D, see below) undergoes further development with the increase of the magnetic field strength forcing out the 
symmetric phase with m — 0. The order parameter oscillation visible in the diagrams when the chemical potential is 
high enough is a phenomenon typical for the model 5ft- Glj, such behavior is generally inherent in cold maiiy-body 
quantum systems in a magnetic field, the fact known since the studies on the de Haas-van Alphen effect [62, Ifij] • 

In general, the NJL model is known to give rise to three distinct phases [6J, [6a] : a symmetric massless phase A 
with no chiral condensate and two chirally broken massive phases B and C, the latter being a phase with a nonzero 
matter density p, whereas p = in phase B (the B ^ C transition occurs for /x > to). Massive phases B and C are 
now spatially nonuniform when H > 0. There is also a new phase with a strong condensate inhomogeneity (retaining 
the presence of DCDW in the i7 — > limit studied in Ref. [20|), we denote it with the symbol D. The position of 
the phases described above in case of if = is illustrated in Fig. [Sfa) and their evolution with the increase of the 
magnetic field strength is shown in Fig. [TJa). The transitions between the phases under consideration are first order 
since the order parameters are discontinuous except for the B ^f C transition when H ~ which is second order 
being a singular point in the diagram [6l[ . It should be noted that the D —>■ C transition occurring in a magnetic field 
strong enough actually belongs to a series of typical order parameter oscillations visible, e.g., in Figs. [2J:, [2Ji. There 
is no significant physical difference between C and D in that region so we consider it as a crossover area plotting the 
corresponding transition with a dotted line, and we only plot a solid line between C and D in the region where these 
phases can be distinguished clearly with a noticeable change in their physical properties [see, e.g.. Figs. [3^, [SJs]. The 
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FIG. 3: The order parameters as functions of the chemical potential for various values of the magnetic field strength at zero 
temperature, and G = 6. All quantities are dimensionless. 



position of the end point separating the solid and the dotted segment is therefore not fixed precisely and should be 
chosen judiciously. In the most general case, one may consider an infinite series of phases {^n}, {Cn}, and {D,i\ 
when a magnetic field is present with phase transitions corresponding to the order parameter oscillations mentioned 
above (like it is done in Refs. |59l - l6l| . see also a recent study in Ref. [33|)- However, since these oscillations are small 
in their relative magnitude and tend to be smeared out with finite temperature taken into account, we consider such 
series as single phases, and in this reasonable approximation, this situation may be treated as a crossover between 
C and D. The main result we have obtained is that there is a nonzero h in all phases when H > except for the 
symmetric one and the case of /i = 0. Smooth and linear growth of b with the increase of H and /i is inherent in 
phase B. Symmetric phase A now occupies a limited area on the diagram. 

We have also examined the case of a subcritical G in addition to the strong-coupling regime. The results for G — 3, 
T = are presented in Figs. [T)d, 21 The magnetic field is known to be a catalyst of the spontaneous chiral symmetry 
breaking both in renormalizable and nonrenormalizable (NJL-like) theories (see, e.g., Refs. [6a - [75| and also Refs. [53- 
|61|). the latter demonstrating the emergence of a dynamic fermion mass for arbitrary small values of the coupling 
constant. This effect is present in our case as well. The phase diagram structure obtained for our model is similar 
to that derived in Ref. [61| for G < Gc and phase B exhibits the same behavior of the order parameter b growth as 
described above. Thus, DCDW formation is preferable for the system in a wide range of the coupling constant. 



VI. CONCLUSIONS 



The calculations performed in the framework of the NJL model have shown that the presence of an external 
magnetic field favors the formation of a spatially nonuniform chiral condensate (in the form of a dual chiral density 
wave) in dense quark matter at low temperatures. This means that there exists a critical magnetic field strength He 



10 



(a) /u = 0.06 



(b) /J = 0.2 



1.0 r 



0.8 - 



0.6 - 



0.4 




1.0 



0.8 



0.6 



0.4 



1.0 



0.8 



0.6 



0.4 



0.2 



0.0^ 
0.8 0.0 




(d) -i/eH = 0.6 



0.2 



0.4 



0.6 



0.8 



FIG. 4: Samples of the order parameter dependence on the external conditions at zero temperature, G = S. All quantities are 
dimensionless. 



such that one of the spatially nonuniform DCDW phases {B, C, or D) emerges in the system when H > He both 
for supercritical and subcritical values of the coupling constant for arbitrary nonzero values of the chemical potential 
/i and T — 0. For example, if G > Gc then it is easy to see that He > for the range of /x corresponding to the 
symmetric phase A, whereas He = for other values of /i [see Fig. [ija)]. On the contrary, in the case of subcritical 
G < Ge, the quantity He is nonvanishing for all /i > [see Fig. [Ifb)]. One can verify that the effect of the chiral 
condensate spatial modulation is mainly due to the particle and antiparticle energy spectra asymmetry induced by 
the presence of DCDW in our model; if one drops out the contribution of the distorted lowest Landau level (LLL) 
to the thermodynamic potential of the system, the phenomenon of the condensate wave vector being nonzero in the 
massive phases of the model except for D will be lost and phase D will be far less stable occupying a small area on 
the phase diagram. 

As discussed in Ref. t20| (see also Ref. [!/£]), linear growth of the condensate wave vector with the increase of the 
chemical potential is generally inherent in one-dimensional systems and this is in agreement with the dimensional 
reduction phenomenon occurring for fermions in a strong magnetic field |7l|-[73| (see also Ref. [76| for the case of 
chromomagnetic fields); this behavior of the order parameters is actually related to the specific properties of the LLL. 
A singular role of the LLL and its impact on physical phenomena in various problems concerning dense matter and 
symmetry breaking is pointed out in a number of studies, see, e.g., Ref. [53[ and also a recent discussion on the chiral 
magnetic effect in Ref. J77J . 

In this paper, we have only reported our results for the case of cold quark matter, and the role of finite temperature 
is to be studied in our forthcoming publications. There are also other interesting subjects left beyond the scope of 
our paper. Since quark matter may possess its own magnetization (see, e.g., Refs. |78l - [8l| and also Refs. [12,13), ^ 
challenging self-consistent problem may arise with the magnetic field being generated dynamically. One should also 
consider the color superconductivity phenomenon possible along with the chiral density waves formation; the results 
obtained in such generalized models seem to be less regularization dependent 22]. A nonzero quark current mass 



11 

should be taken into account as well. Besides, concerning the ground state spatial configuration of the NJL model, 
it has been argued that domain walls may be more preferable than chiral density waves at least in the absence of 
external gauge fields [2^. At the same time, as it may be concluded from our calculations, one would expect that a 
strong magnetic field favors the formation of DCDW. Thus, there should exist some solution interpolating between 
the two extremes in the intermediate region of the magnetic field strength, possibly being similar to the solution 
discussed in Ref. [82] . On the other handja competing mechanism for domain walls formation in a strong magnetic 
field has also been discussed in literature [42| . Therefore, the problem of the preferred ground state spatial structure 
requires further theoretical investigation but, in general, it has been shown that an external magnetic field induces 
condensate inhomogeneity, in the form of DCDW or some more preferable configuration. Another subject of research 
is obtaining analytical expressions for the order parameters as functions of the external conditions in a weak magnetic 
field at least in some special cases using an approach similar to that adopted, e.g., in Ref. |83| . 

Our concluding remark is that real existence of a spatially nonuniform chiral condensate in nature is yet an open 
question since theoretical results related to the problem are generally model and approximation dependent. Unfor- 
tunately, exact QCD calculation of its production is impossible, since this is an infrared phenomenon. Nonetheless, 
we believe that the theoretical research of this kind of nonperturbative effects will yield our better understanding the 
properties of strongly interacting matter. 



Appendix: Regularization of Q'^ at the n — Level 



Let us consider the contribution of the n = energy level in expression (j2ip : 

+ 00 



K"'= I dpY,{\E-lJ-\-\E\)e{K' -\E\), E = e^m^+p^ + b, 





(A.l) 



where we have utilized the parity of E with respect to p under the integral and omitted the constant factor — ei?/(27r)^. 
One can easily prove the following formulas being valid for sufficiently large values of P: 



I 



dp 



\/ vn? +p^ + a 



I{Q,P) + Pa, a>-|TO|, 

-I{0,po) + lipo, P) - 2poa + Pa, a < -\m\, 



dp 



-\/ m^ + p^ + a 



I{0,P)-Pa, a<|m|, 

~I{0,Po)+I{po,P) + 2poa~Pa, a> |m|. 



where po — yoF^-rr? and 



I{pi,P2) = dpy/ 



Pi 



rri^ , I P2 + \J"m?- + p\ 



2 ^ 2 V 2 \ pi + ^"^ +Pi 



Let us apply the above formulas to evaluate K'^^ . We have to choose different integration limits for different values of e 
due to the energy cutoff and the spectrum asymmetry: P\ — •\/(A' — 6)^ — m? for e = +1 and P^ = ^(A' + 6)^ — m^ 
for e — —1. Consider the following expression: 



dp 



\/ rrp- + p2 _j_ a 



F2 



dp 



- \/ni^ + p^ + I 



/(po. Pi) + /(po, Pi) + (Pi - Pi)a + 2poa, \m\ < a, 
/(0,Pi) + /(0,P2) + (Pi-P2)a, -\m\ <a<\m\ 

I{po,Pi) + Iipo,P2) + (Pi - P2)a - 2poa, a < -\m\. 



(A.2) 
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Let us compare this result to one obtained with the help of a trivial momentum cutoff (i.e., without the factor 
9{K' — \E\) but with a common upper limit in the integrals). Assuming Pi = P2 = A', we get 



A' 



J" 



dp 



\/rn? +\ 



■ p^ + a 




dp 



- \lm?- + p2 + a 



2/(po,A')+2poa, |to| < a, 

2/(0, A'), -|m|<a<|m|, 

2I{pQ,A')~2poa, a<-\m\. 

The difference of the above quantities is 
Thus we have 



(A.3) 



L^cn Tcn Tcn rn 

A —J \a^b-Li — -/ |a^6 — -/ 



|a— 6— /i 



-J" 



|a=6 + KP2 -Pl)= K'^°'^ + ^i{P2 - Pi), 



where we have introduced the symbol X™°™ to denote expression (jA.l[) regularized with a trivial momentum cutoff 
instead of the 6'(A' - |£'|) factor: 



A' 



K^°^ = j dpY,{\E-n\-\E\ 



(A.4) 



Since 



we finally get 



lim (F2 -Pi) = 26, 

A'— J- 00 



+00 A 

j dpY^{\E-^,\-\E\)e{K'-\E\) =JdpY,i\E-^,\-\E 



A'^-CXD Q 



2^6. 



(A.5) 



A'— >oc 



We should note that expression (jA.4[) for iif™°™ is convergent and well-defined in the limit A' ^ cx). Taking sufficiently 
large values of p, one has 



y I ley'm'^ + p'^ + b — ^ — €\Jm'- +p2 + 6 j = 
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